Experimental observation of the Bogoliubov transformation for a Bose-Einstein 

condensed gas 
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Phonons with wavevector q/h were optically imprinted into a Bose-Einstein condensate. Their 
momentum distribution was analyzed using Bragg spectroscopy with a high momentum transfer. 
The wavefunction of the phonons was shown to be a superposition of +q and — q free particle 
momentum states, in agreement with the Bogoliubov quasiparticle picture. 
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The pioneering work of Bogoliubov in 1947 constitutes 
the first microscopic theory that attributes superfluidity 
to Bose-Einstein condensation jjj. As described by Ein- 
stein noninteracting bosons condense Q , but they are not 
superfluid. However, Bogoliubov showed that with weak 
interactions the condensate will exhibit superfluidity. Re- 
pulsive interactions change the elementary excitations at 
long wavelengths from free particles into phonons which, 
according to the Landau criterion, leads to superfluidity 
The main step in the non-perturbative treatment is 
the Bogoliubov transformation 
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which expresses the creation and annihilation operators 
b q and b q for Bogoliubov quasiparticles in terms of cre- 
ation and annihilation operators a q and a q for free par- 
ticles in momentum states q. This transformation also 
plays a crucial role in general relativity, where it connects 
the particle operators in different reference frames Q . In 
this paper we experimentally verify the Bogoliubov trans- 
formation by generating quasiparticles with momentum 
+q and observing that they are a superposition of +q 
and — q momentum states of free particles. 

Following a recent suggestion of Brunello et al. ||, we 
perform an experiment that combines the two regimes 
of Bragg spectroscopy, where the momentum imparted 
to the atoms is either smaller || or larger |?|] than the 
sound velocity c times the mass m of the atoms. An op- 
tical lattice moving through the condensate at the sound 
velocity imprinted phonons with wavelengths equal to 
the spatial period of the lattice Q. This was accom- 
plished by intersecting two laser beams at a small angle 
and choosing their frequency difference to be equal to 
the phonon frequency. The momentum analysis of the 
phonon wavefunction was performed by a second Bragg 
pulse consisting of counterpropagating laser beams that 
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transferred a large momentum Q (two photon recoils) to 
atoms with initial momentum p. In the limit Q 3> p, 
mc, the resonance frequency v is equal to the kinetic en- 
ergy transferred: hv — Q 2 /2m + Qp/m, where the second 
term is simply the Doppler shift [ft |], 9. The resulting 
frequency spectrum shows three peaks corresponding to 
the three momentum components (0, +q and —q) of a 
condensate with phonons (see Fig. |2[) M . 

The Bogoliubov spectrum for the energy of elementary 
excitations is ifTol 



£(?) = 



q 2 c 2 



2m 



(2) 



where c = V4irh 2 an/m, a is the scattering length, and 
n is the density. The amplitudes u q and v q in Eq. |lj are 
given by 



= ± 
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The non-trivial aspect of the Bogoliubov transforma- 
tion manifests itself only at low momenta (q <C mc), 
where both amplitudes u q and v q are significant. In 
this regime, the excitations are characterized by u ~ 
— v ~ y / mc/2q 3> 1 and e(q) ~ qc. Such excitations 
are phonons, each involving many particles moving in 



both directions. At high momenta (q 3> mc), 
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v q ~ and e(q) = q /2m + mc . These excitations are 
free particles with an energy shift equal to the chemical 
potential [i — mc 2 . 

The phonon creation operator b +q (Eq. |l|) is a super- 
position of the creation operator a} +q for particles moving 
in the +q direction and the annihilation operator a_ f/ for 
particles moving in the —q direction, yet the creation of 
a phonon in the condensate implies an increase of par- 
ticles moving in both the +q and — q direction. Indeed, 
simple operator algebra shows that a condensate with I 
excitations b+ q / VV-l^o) contains lu q + v 2 free particles 
moving with momentum +q and (/ + l)v 2 free particles 
with momentum —q. In its ground state the conden- 
sate contains v 2 pairs of atoms with momenta +q and 



—q. These pairs constitute the quantum depletion in the 
condensate wavefunction [[10| 
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are in the q = mo- 
mentum state. When a\_ q and a_ g act on 'J'oj 
terms with large occupation numbers j are enhanced: 
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addition, — q and +q atoms appear only in pairs. 



In 
To- 
gether, these two effects cause both a) +q and a_ g to in- 
crease the number of atoms in both the — q and +q states. 

The experiments were performed with condensates of 
3 x 10 7 sodium atoms in a magnetic trap with radial 
and axial trapping frequencies of 37 Hz and 7 Hz, re- 
spectively jllj. The condensate had a peak density of 
1.0 x 10 14 cm -3 corresponding to a chemical potential of 
fj, = h-1.5 kHz, a sound velocity of c = 5 mm/s, and a 
Thomas-Fermi radial radius of 32 /im. The Bragg beams 
for the optical lattices were generated from a common 
source of laser light 1.7 GHz red-detuned from the 3Si/ 2 
\F = 1) to 3P 3/2 \F' = 0, 1,2) transitions. The lattices 
were moved radially through the cloud. This was done 
to avoid the high collisional density along the axial di- 
rection, where outcoupled atoms would have undergone 
elastic collisions with the condensate ]l^| . The Bragg 
beams for imprinting the phonons were at an angle of 32 
mrad with respect to one another, resulting in a lattice 
spacing of ~ 9 /im and a recoil momentum of q — to- 1.9 
mm/s. The beams propagated at an angle of 0.4 rad 
with respect to the longitudinal axis of the condensate 
and were linearly polarized perpendicular to it. The fre- 
quency difference ("excitation frequency") between the 
two beams was chosen to be 400 Hz, corresponding to 
the frequency of the phonons. The excitation pulse had 
to be long enough to ensure sufficient frequency resolu- 
tion in order to selectively excite +q phonons, and no 
— q phonons (see below). However, the excitation pulse 
also had to be shorter than the transit time of phonons 
through the condensate, since phonons accelerate when 
they move through regions of varying density. We chose 
small-angle Bragg beams 3 ms in duration with inten- 
sity 0.05 mW/cm 2 , corresponding to a two-photon Rabi 
frequency of 50 Hz. 

The momentum analysis of the phonons was performed 
with two counterpropagating beams that imparted a re- 
coil momentum of Q = m-59 mm/s onto the outcoupled 
atoms. The beams were polarized parallel to the longi- 
tudinal axis of the condensate to suppress superradiant 
A frequency difference ( "probe frequency" ) 
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of 100 kHz between the two beams corresponded to the 
kinetic energy needed for atoms initially at rest to reach 
this recoil momentum. To gain recoil momentum +Q, 
atoms with initial momentum +q were resonant at 107 
kHz (Fig. |l|a). Atoms with momentum — q were resonant 
at 94 kHz (Fig. [Ob). In our experiment, a retroreflected 
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FIG. 1: Momentum distribution of a condensate with 
phonons. After imprinting +q phonons into the condensate, 
momentum analysis via Bragg spectroscopy transfers a mo- 
mentum ±Q (two photon recoil) to the atoms. Absorption 
images after 40 ms time of flight in (a), (b), and (c) show the 
condensate in the center and outcoupled atoms to the right 
and left for probe frequencies of 94, 100, and 107 kHz, respec- 
tively. The small clouds centered at +q are phonons that were 
converted to free particles. The size of the images is 25x2.2 
mm. (d) The outlined region in (a) - (c) is magnified, and 
clearly shows outcoupled atoms with momenta Q ± g, imply- 
ing that phonons with wavevector q/h have both +q and — q 
free particle momentum components. 



beam containing both optical frequencies resulted in two 
optical lattices moving at the same speed in opposite 
directions. This led to simultaneous outcoupling of +q 
and — q atoms in opposite directions. These large-angle 
Bragg beams were pulsed on for 0.5 ms. The probe pulse 
had to be long enough to selectively excite atoms with 
±q momentum, but also shorter than h/fi for the mean 
field energy to be negligible during the readout. Each 
frequency component had an intensity of 0.7 mW/cm 2 , 
corresponding to a two-photon Rabi frequency of 700 Hz. 
Subsequently, the trap was turned off and a resonant ab- 
sorption image was taken after 40 ms of ballistic expan- 
sion. 

Fig. 1 shows typical absorption images for various 
probe frequencies. The quasiparticle nature of the 
phonons was directly evident (see outlined region) in the 
time-of-flight distribution through the presence of the 
peaks at momenta ±q + Q. These peaks had well de- 
fined momentum because the outcoupled atoms left the 
condensate quickly (during this time the atoms' velocity 
changed by less than the speed of sound Q ) . This "pho- 
tograph" of the Bogoliubov transformation is the central 
result of this paper. 

We now discuss the different momentum components 
distinguishable in Fig. 1: 

(1) The original condensate is in the center. 

(2) The condensate is asymmetrically extended to- 
wards positive momenta. This is due to imprinted 
phonons of momentum q that were converted to free par- 
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FIG. 2: Bragg spectrum of a condensate with (•) and without 
(o) phonon excitation. The number of outcoupled atoms vs. 
the probe frequency is shown. The excitation of phonons cre- 
ates sidebands on both sides of the central condensate peak. 
Around 100 kHz, the absorption images were saturated due 
to high optical density. 



tides during the ballistic expansion. Previously this sig- 
nature was used to determine the structure factor of a 
condensate Q. There is no momentum component in 
the opposite direction, because we did not create —q 
phonons (see below) . The smearing of the observed mo- 
mentum distribution may be caused by acceleration due 
to the inhomogeneous density distribution of the conden- 
sate when the trap was switched off. 

(3) The two components with momentum +Q and — Q 
are atoms outcoupled from the condensate at rest. The 
symmetry of the ±Q peaks, as well as the position of the 
main condensate, served as an indicator as to whether the 
condensate was undergoing dipole oscillation during the 
experiment. Images that showed such "sloshing" (caused 
by technical noise) were excluded from further analysis. 

(4) Atoms at +q + Q (Fig. lc) and +q - Q (Fig. la) 
are coupled out from the +q component of the phonons. 

(5) Atoms at -q - Q (Fig. lc) and -q + Q (Fig. la) 
are coupled out from the — q component of the phonons. 

Quantitative information was obtained by scanning the 
probe frequency and measuring the number of outcou- 
pled atoms in the three peaks around — Q (See Fig. 1). 
Without phonons only the condensate peak is observed. 
The excitation of phonons at wavevector +q creates mo- 
mentum sidebands at momentum ±q. The — q peak is 
expected to be smaller by a factor Vg/iig. This effect is 
evident in Fig. 2, but with a poor signal-to-noise ratio. 

The excitation of the phonons was characterized by 
scanning the frequency difference of the small angle 
Bragg beams (Fig. 3a) and measuring the number of out- 
coupled atoms at +q — Q. The probe frequency was kept 
at 94 kHz because the momenta of the excited phonons 
remained fixed. The observation of two distinct peaks 
(corresponding to +q at 400 Hz and — q phonons at - 
400 Hz, respectively) confirms that there was sufficient 
resolution to excite only +q phonons and suppress the 
off- resonant excitation of — q phonons. The higher peak 
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FIG. 3: Phonon excitation spectrum. Atoms with initial mo- 
mentum +q were detected by setting the probe frequency to 
94 kHz and measuring the number of atoms with momentum 
— Q + q. The two peaks reflect that phonons with +q/h and 
—q/h wavevectors have free particle components with momen- 
tum +q. Spectra were taken at: (a) high density (l.OxlO 14 
cm -3 ) and (b) low density (0.5xl0 14 cm -3 ). The solid line, a 
fit to the sum of two Gaussians, is intended to guide the eye. 
The dashed line is the theoretical prediction. 



represents the u-component of +q phonons excited at 
positive excitation frequencies. When the excitation fre- 
quency became negative, — q phonons were excited. This 
second peak represents the w-component of the phonons. 

According to Eq. ||| the ratio of the two peaks, v^/Ug, 
should be smaller at lower density. This is confirmed 
in Fig. 3b, which shows the excitation spectrum at low 
density. The density was lowered by a factor of 2 by 
weakening the axial trap frequency to 4 Hz and reducing 
the number of atoms in the condensate by a factor of 3. 
The scatter in the data could be due to residual sloshing, 
shot-to-shot fluctuations in the size of the condensate, or 
non-linear effects (see below). 

The dashed line is a the theoretical prediction with 
no free parameters (except the vertical scale). It was 
obtained by integrating Eq. [}] over the inhomogeneous 
density distribution of the condensate and by accounting 
for the finite length of the square excitation pulse, which 
broadens the line and causes the extra sidelobes around 
900 Hz. The theory assumes the validity of a perturbative 
approach, i.e. that the excitation pulse is weak. 

Ideally, both Bragg pulses should affect only a few per- 



4 



cent of the atoms, as in previous experiments l| Q . How- 
ever, because our signal was the product of the two out- 
coupling efficiencies, we needed to work at much higher 
outcoupled fractions. We estimate that the excitation 
pulse transferred 10 of the condensate atoms into the 
phonon state, and that the probe pulse outcoupled 40% 
of these atoms on resonance. 

Both Bragg processes (stimulated Rayleigh scattering) 
should depend only on the product of the two intensities. 
Therefore, changing the sign of the excitation frequency 
should not affect the number of phonons generated (as- 
suming a stationary condensate). However, when the 
sign of the excitation frequency was changed in our ex- 
periment, the observed number of phonons differed. This 
is most likely due to superradiance , which is sensitive 
to the individual intensities. Therefore, the asymmetry 
in phonon number was eliminated by ensuring equal in- 
tensities in the two beams. Additionally, there was a 
substantial loss (^50%) of condensate atoms due to su- 
perradiant Rayleigh scattering. Both superradiant effects 
could be further suppressed using light further detuned 
from resonance. However, in our case this would have 
required an additional laser. Given these experimental 
limitations, the agreement between experiment and the- 
ory in Fig. 3 is satisfactory. 

In this work we have used large-angle Bragg pulses to 
analyze the momentum structure of the phonon wave- 



function. In principle, this could have also been achieved 
by removing the mean-field interaction within a time 
h/e{q) and then probing the velocity distribution of the 
particles. This is not possible with ordinary ballistic ex- 
pansion because the reduction of the interactions is too 
slow, taking place on a time scale of the trapping period. 
However, the use of a Feshbach resonance [|l5| would pro- 
vide an effective method for suddenly reducing the mean- 
field interaction. 

In conclusion, we have experimentally analyzed the 
phonon wavefunction in a Bose-Einstein condensate. Fol- 
lowing recent theoretical work the two-component 
character of Bogoliubov quasi-particles was observed in 
the frequency domain (Fig. 2). In addition, the momen- 
tum components of the phonon wavefunction were dis- 
criminated by their final momenta after the probe pulse 
(Fig. 1). By combining momentum and frequency selec- 
tivity, we were able to directly "photograph" the Bogoli- 
ubov transformation (Fig. Id) , demonstrating the power 
of Bragg spectroscopy to analyze non-trivial wavefunc- 
tions. This method may be also applicable to study- 
ing the many-body and vortex states Q of dilute atomic 
Bose-Einstein condensates. 

This work was funded by ONR, NSF, ARO, NASA, 
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